Non-equilibrium Enhancement of Cooper Pairing in Cold Fermion Systems 
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Non-equilibrium stimulation of superfluidity in trapped Fermi gases is discussed by analogy to the 
work of Eliashberg [G. M. Eliashberg, in "Nonequilibrium Superconductivity," edited by D. N. Lan- 
genberg and A. I. Larkin (North-Holland, New York, 1986)] on the microwave enhancement of su- 
perconductivity. Optical excitation of the fermions balanced by heat loss due to thermal contact 
with a boson bath and/or evaporative cooling enables stationary non-equilibrium states to exist. 
Such a state manifests as a shift of the quasiparticle spectrum to higher energies and this effectively 
raises the pairing transition temperature. As an illustration, we calculate the effective enhancement 
of Cooper pairing and superfluidity in both the normal and superfluid phases for a mixture of 87 Rb 
and 6 Li in the limit of small departure from equilibrium. It is argued that in experiment the desir- 
able effect is not limited to such small perturbations and the effective enhancement of the pairing 
temperature may be quite large. 
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I. INTRODUCTION 

The difficulty of observing quantum coherent phases in 
cold gases highlights the need to overcome low transition 
temperatures. In addition to Bose-Einstein condensa- 
tion, there has also been recent interest in the genera- 
tion of Fermi superfluids through the BCS pairing mech- 
anism 1]. This phenomenon is more difficult to observe 
due to prohibitively low transition temperatures @, Q 
though the problem may be partially surmounted by use 
of Feshbach resonances [H, 0, [5| . Non-equilibrium effects 
can also be used to control and effectively cool such sys- 
tems. However, this is an unexplored area of research by 
comparison. 

Since the 1960's, it was known that superconductivity 
could be stimulated by radiation in microbridges [6[ . In 
1970, Eliashberg explained this effect as an amplification 
of the gap parameter by means of a stationary nonequi- 
librium shift in the quasiparticle spectrum to higher en- 
ergies brought on by the radiation 0, H[. Over the 
next decade, his theory found experimental acceptance 
through the enhancement of critical currents and temper- 
atures in Josephson junctions [9| and thin films [l(| • At 
the same time, other nonequilibrium stimulation meth- 
ods were developed [ll| with more recent reports of en- 
hancements of the superconducting critical temperature 
by up to several times its equilibrium value [12| . Il3j . With 
the present interest in the application of the BCS model 
of superconductivity to trapped atomic Fermi gases 
Il5j , nonequilibrium effects represent an attractive way to 
magnify the quantum properties of these types of super- 
fluids. 

As in superconductors, the BCS order parameter Ao 
for cold fermionic gases obeys a self-consistency equation 
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where £k = £k — £f is the quasiparticle dispersion cen- 
tered on the Fermi energy e-p, Ao is the BCS gap, and 
rtk is the quasiparticle distribution function. The con- 
stant A has the form A = —Ana^i/mF where a-j-j, is the 
negative s-wave scattering length for collisions between 
hyperfine states and mp is the mass. At equilibrium, 
n k * s ^ ne Fermi-Dirac distribution function, and the 
only way to increase Ao is either to increase the inter- 
action strength or to lower the temperature. However, 
there exists a wide class of stationary nonequilibrium dis- 
tributions, rik, such that Eq. ([T]) is still valid and has 
solutions with enhanced order parameters. Indeed, if 
a quasistationary non-equilibrium distribution is created 
that is different from the canonical Fermi-Dirac function, 
<5rik = «k — n ^ D \ then according to the weak-coupling 
BCS Eq. ([1]), it effectively renormalizes the pairing inter- 
action and transition temperature as follows: 
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where here and below 
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the density of states at the Fermi level, Tc 
Ef exp {— 1/ [i/(ep)A]} is the weak-coupling BCS transi- 
tion temperature in equilibrium, and we also introduced 
the dimensionless parameter x — ~Sk^ n k/ [i / (ef)Ej s \. 
For many non-equilibrium distributions, % > 0, and this 
yields an effective enhancement of T c and/or A. We note 
that even though our theory below and that of Eliash- 
berg are limited to small deviations from equilibrium, 
with \x\ "C 1, this does not imply a limitation in experi- 
ment, where this parameter can be large. For such large 
deviations from equilibrium, the weak-coupling BCS ap- 
proach and Eqs. |T]) and {2J may not be quantitatively 
applicable, but the tendency to enhance pairing may re- 
main. Therefore the proposed underlying mechanism 
may lead to substantial enhancement of fermion pair- 
ing and superfluidity. We also emphasize here that cold 
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FIG. 1: Bragg Potential: A moving lattice with wave vector 
q = k2 — ki can be formed in the region of the Bose-Fermi mix- 
ture through the interference of two lasers with differing wave vec- 
tors and frequencies. By adjusting the parameters of this non- 
equilibrium perturbation, one can achieve states with enhanced 
superfluidity. 



atom systems offer more control in creating and manipu- 
lating non-equilibrium many-body quantum states than 
that available in solids. Specifically, we will show that 
while it was impossible to drive a metal from the normal 
to the superconducting phase by irradiation, it is indeed 
possible to drive the equivalent transition in cold gas by 
utilizing this additional control. 

In this paper, we propose a theory of nonequihbrium 
stimulation of fermion pairing by considering the effect of 
Bragg pulses [ijllIII as shown schematically in Fig. [T] on 
a harmonically trapped gas of fermions in the Thomas- 
Fermi approximation [20(. The heating induced by the 
external perturbation is dumped into an isothermal bath 
of trapped bosons via collisions, but this is not necessary 
in general. The pairing enhancement is calculated for a 
typical mixture of 87 Rb and 6 Li. It depends on the state 
of the gas at equilibrium. In the superfluid phase, Eliash- 
berg's requirement that the frequency of the perturbation 
be less than twice the equilibrium gap (Tiw < 2A ) en- 
sures that the pulse does not effectively heat the system 
by producing more quasiparticles with energies e ~ Aq. 
Though this requirement cannot be satisfied in the nor- 
mal phase, where Ao = 0, the independent tunability 
of both the momentum and energy of the Bragg pulse 
allows us to protect the system from effective heating 
through energy conservation. This avenue, which was 
not available in the context of superconductors, effec- 
tively provides a means to "sharpen" the Fermi step (or 
even create a discontinuity at a different momentum), 
thereby enhancing fermion pairing. 



II. MODEL 

In our problem, we assume optically trapped bosons in 
thermal contact with fermions that occupy two hyperfine 
states |t) and |4-}. This system has a Hamiltonian of the 
form H = Ho + Hi where the noninteracting part Hq is 



given by 
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where for brevity, we introduced the subscript p = 
B, Ff, F± which labels bosons (p — B) and fermions in the 
two hyperfine states: "up" (p = Ff) or "down" (p = Fj,) 
and mF t = wf 4 = tof is the fermion mass and ms is 
the mass of the bosons. We assume that fermions in 
either hyperfine state feel the same trapping potential. 
Thus, V p (r) is given by Vf,b (r) = ^ m F,B&F,B r2 where 
the subscript F (B) refers to fermions (bosons). There is 
also an interaction Hamiltonian Hi which has the form 
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with g piP2 being the strength for s-wave collisions be- 
tween the particles labelled by pi,P2 — {B,Ff,F±}. 
While Pauli exclusion requires that <?F t F t = 9f 1 f 1 = 0, 
an attractive coupling gF t F_ l = 9f±f^ < will lead to 
BCS pairing. A nonzero interaction gpB between bosons 
and fermions is required for thermalization between the 
two populations. We need put no other restrictions on 
9piP2 i but our desired effect will be easier to observe ex- 
perimentally with some other constraints. For instance, 
requiring that Qwb < will raise the BCS condensation 
temperature [2 11 ] while a larger qbb > facilitates ther- 
malization between bosons and fermions. 

To proceed further we use the gap equation, Eq. |T]), 

where we have — n^ B ^ at equilibrium. In the Thomas 

Fermi approximation, the transition temperature, Tc°^ 
given by @ 
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where kp is the Fermi wave-vector and 7 ~ 0.577 ... is 
Euler's constant, the scattering length a-fj. is a simple 
combination of the coupling strengths g PlP2 ■ 



A. Nonequilibrium Enhancement 

It is possible to create distributions that lead to larger 
order parameter and effective condensation temperature 
by weakly perturbing the trapped fermions. Specifi- 
cally, we affect a Bragg pulse (Fig. [TJ by illuminating 
the fermions with two lasers, which are both largely de- 
tuned from any fermionic transition. In what follows, we 
assume the lasers to be even further detuned from any 
bosonic transition such that we may ignore the effect of 
the Bragg pulse on the bosons. The interaction of the 
fermions with these lasers is described by the addition of 
a term 
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to the Hamiltonian where q and ui represent the dif- 
ference in wavevectors and frequencies between the two 
lasers [l^. Now, following Schmid and the general 
argument in the introduction [see, Eq. ([5])]. we introduce 
the function, Sn^, which describes departure from equi- 
librium. While Eliashberg assumed in |8| that the impu- 
rity concentration was high enough in metals such that 
momentum relaxation happened at a much faster rate 
than energy relaxation, we shall not make this assump- 
tion. As such, Srik need not be isotropic although this 
requirement is easily included in our model. The cor- 
responding term, x/^(ef) from Eq. is added to the 
right side of the gap equation ^ and leads to a new solu- 
tion, A > Ao, for the order parameter at the same tem- 
perature (note that in the non-equilibrium situation, the 
notion of a temperature of the Fermi system is undefined, 
and here by temperature we imply the original temper- 
ature and that of the Bose bath). For T - T c (0) « T c (0) , 
Eq. |T]) can be cast in the form of a Ginzburg-Landau 
equation, which, including the nonequilibrium term, be- 
comes 

(l»!l-*)A + t (i§) 2 4 = „, (4) 

where we assume that the coefficient in the cubic term of 
Eq. ((31) is only weakly affected by the perturbation and 
use its standard BCS value [C(z) is the Riemann zeta- 
function] b = 7C(3)/(8tt 2 ) » 0.107. . . (see also, Ref. H). 
Eq. (j4j nominally leads to an exponential enhancement 

of the effective critical temperature: T^ ff = e-ri ' , if \ 
is positive. 

B. Kinetic Equation 

To calculate this enhancement for the Bose-Fermi mix- 
ture, we shall balance the Boltzmann equation for in 
the spirit of Eliashberg, including both a collisive contri- 
bution and that from Bragg scattering. 

flu = Icoll [™k] + ^Bragg [™k] ■ (5) 

For small departures from equilibrium, we can lin- 
earize the collision integral in Sn^ and use the 1/r- 
approximation: 7 co n [n] — —Sn/T , where To is the quasi- 
particlc lifetime. In our system, this lifetime will be dom- 
inated by the inelastic collision time between bosons and 
fermions. We may estimate this time as in 0, HH by 
means of a 1/tq = nav approximation where n is the 
boson density at the center of the trap, a is the constant 
low temperature cross-section for boson-fermion scatter- 
ing, and v is the average relative velocity associated with 
the collisions between bosons and fermions. 

Note that there exist other contributions to the col- 
lision integral, in particular those coming from fermion 
collisions. Our model assumes point-like interactions be- 
tween fermions: Such interactions can be separated into 
interactions in the reduced BCS channel, which involve 



particles with opposite momenta that eventually form 
Cooper pairs and other types of scattering events, which 
give rise to Fermi liquid renormalizations on the high- 
temperature side and superconducting fluctuations on 
the BCS side. Note that dropping off the latter terms 
would lead to an integrable (Richardson) model that 
does not have any thermalization processes and there- 
fore the collision integral for its quasiparticles must van- 
ish. In thermodynamic limit this model is described by 
BCS mean-field theory perfectly well and so we can say 
that the pairing part of fermion interactions is already 
incorporated in our theory. Of course, including fermion- 
boson collision and the second type of fermion interac- 
tion processes break integrability and lead to two types 
of effects: First, such interactions lead to Fermi liquid 
renormalizations of the effective mass and the quasipar- 
ticle ^-factor. However, these effects are not germane 
to the physics of interest, and we may assume that all 
relevant corrections arc already included and treat our 
system as that consisting of Fermi liquid quasiparticles. 
However, there is of a course a second dissipative part 
coming from interactions, such as those due to bosons al- 
ready included into tq and quasiparticle scatterings and 
decay processes due to non-BCS fermion-fermion colli- 
sions. Similarly to the work of Eliashberg, we will as- 
sume that the latter contribution to the collision integral 
is less significant than Tq 1 due to the bosons. Fermi 
Liquid quasiparticles are exactly defined precisely on the 
Fermi sphere, but they have a finite lifetime due to de- 
cay processes elsewhere. By not including the lifetime 
Tk of a quasiparticle at momentum k in our linearization 
of the Boltzmann equation, we have implicitly assumed 
that Tk ^> t . Because — oc (Trk B T) 2 + (e k — Ep) 2 in a 
Fermi Liquid (26l |. there will always be an energy region 
where this assumption will indeed be true for low temper- 
atures. The most important contribution to the integral 
in the expression for \ comes from states for which £k is 
within ksT of ep. As such, if we require that T <C Tp 
and hu <C ep, then our linearization of the Boltzmann 
equation with respect to To will be legitimate for the cal- 
culation of an enhancement of superfluidity. Again, we 
stress the importance of recognizing that the aforemen- 
tioned requirements are necessary only for quantitative 
accuracy of our model. As with Eliashberg's enhance- 
ment of superconductivity, we expect our effect to be 
observable far outside the constrained parameter space 
that is necessary for strict validity of our simple model, 
which provides a proof of principle for using nonequilib- 
rium perturbations to enhance fermion pairing in cold 
atom systems. 

With these caveats in mind, we shall tune the fre- 
quency of our Bragg pulse such that lot$ 1. This 
will ensure that any non-stationary part of the distri- 
bution function will be small Q. Equivalently, we may 
think of this requirement as the statement that the Bragg 
pulse pumps the system out of equilibrium must faster 
than the system relaxes. We may note here that the 
aforementioned assumption that Sn is small also im- 
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plies that x <C 1, thereby diminishing the desired effect, 
Tf = (1 + x)T c (0) - T c (0) . Again, this approximation 
greatly simplifies our theoretical problem by allowing us 
to expand the Boltzmann equation, but it is only a math- 
ematical convenience that represents no impediment to 
an experimentalist looking for striking enhancements of 
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Eq. ([5]) can now be solved for Sn to yield 
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which shows that a stationary nonequilibrium state is 
formed in a characteristic time r . The Bragg part in 
Eq. ([5]), ^Bragg [n], now depends only on n| c FD ' ) and can 
be computed with Fermi's golden rule. When the wave- 
length of the Bragg pulse is much larger than the De- 
Broglie wavelength of the fermions and the reciprocal 
frequency of the pulse is much smaller than the relax- 
ation time (Xp |q| <C 1 and wr >• 1), Fermi's golden 
rule yields 
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The determination of ^Bragg [ n k] aii ows us to find x 

and ultimately Tc° S K The optimal parameters depend 
on whether the fermionic gas is in the superfluid or nor- 
mal phase at the time the Bragg pulse is applied. In the 
former case, the energy gap in the quasiparticle density 
of states protects a Bragg pulse with hcu < 2Ao from 
producing new quasiparticles that will hinder superflu- 
idity. However, in the normal phase, the energy conser- 
vation requirement that £k+q — £k = ihcu and an inde- 
pendent control of q and lo allow us to engineer a pulse 
that ensures that only "thermal" quasiparticles with en- 
ergies e > £f are pushed to even higher energies, while 
the fermions below ep are not affected. 

Thus far, we have assumed that heat is being dissi- 
pated from the fermions into the bosons through colli- 
sions. Because Bose-Einstein condensation inhibits col- 
lisions with fermions and severely reduces thermaliza- 
tion between the two populations, our simple analysis 
depends on the bosons being at a constant temperature 
T greater than their Bose-Einstein condensation temper- 
ature Tbec - 0.94hn B N B /3 [13. Condensation can be 
prevented at temperatures close to the BCS transition 
temperature by having Qb <C &f and tub 3> rap. Treat- 
ing the bosons classically, we expect their temperature to 
increase no faster than ^ = -^QbgUj, which is the energy 
pumping rate due to the Bragg pulse for a specific heat 
C. For a harmonically confined classical gas, we use a 
specific heat given by C = SksNs- If tsragg is the time 
over which the Bragg pulse is turned on, then so long 
as tBragg^ is much less than the temperature of the 
bosons, we may consider the bosonic population to be 
an isothermal bath. This may be accomplished by hav- 
ing a large number of bosons at low density. One may 
be able to avoid the assumption of a bosonic population 
altogether when driving the transition from normal to su- 
perfluid at temperatures above Tj ' by allowing energetic 
particles to leave the trap as in evaporative cooling. We 
shall show later that this is possible because the Bragg 
pulse may be tuned to couple only to particles with en- 



ergies above a threshold energy depending on uj and q. 
For concreteness however, we shall keep the bosonic pop- 
ulation throughout the following section. 



III. NUMERICAL RESULTS 

As an example, we calculate the nonequilibrium en- 
hancement of T c for a trapped mixture of 87 Rb and 6 Li 
under the aforementioned assumptions with 5n <C 1. We 
assume a cloud of 10 5 Lithium atoms and 10 7 Rubidium 
atoms in traps of frequencies rip = 200£ls = 3 kHz cor- 
respondingly. We use scattering lengths 0,3 b = 109a , 
apF = — 2160ao, apB = — lOOao [Il[ and a typical colli- 
sion time To ~ 136 ns as estimated via the 1/tq — nav 
approximation from above. With these parameters, the 
equilibrium BCS and BEC condensation temperatures 
are T c (0) « 0.15 T F = 291 nK and T B ec ~ 23.2 nK. The 
quasiparticle lifetime Tk for the normal fluid is estimated 
as 7k ~ 3 /is for |£k — £f| — via the methods in Refs. 

MM. 



A. Superfluid at Equilibrium 

Let us assume that the system is initially in the su- 
perfluid phase at equilibrium with a low enough tem- 
perature T < T c (0) such that 2A > Twj = 0.15e F - Us- 
ing Eq. (O, we may calculate the stationary distribution 
function (Fig. [5]) for Bragg parameters = 70S! f and 
|q| = 0.1 hp- Some comments on the form of Sn^ are 



necessary. As shown in Fig. 2(a) the first order approx 



imation to rik has unphysical singularities at E — A and 
E = A + to. These are due to first-order transitions to 
E = A+oj from E = A where the quasiparticle density of 
states diverges for a superfluid at equilibrium. The exact 
distribution function, schematically drawn in Fig. 2(b) [ 
has no infinities. Higher orders in the expansion of the 
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FIG. 2: | (a) | The first order approximation to the quasiparticle 
occupation as a function of E = \J S, 2 + Aq for parameters Q(, s = 
70 tip, |q| = O.lkp , and hui = 0.l5ep. The unphysical singularities 
at E = Ao and E = Ao + hui are not included in the calculation of 
X- See text for details. The equilibrium values for this system are 
T = 0.14T F and A = 0.08e F - [(b)]The exact distribution function 
schematically drawn for the same parameters with Ao and hui in 
units of ep. The thick dashed lines represent the occupation at 
equilibrium. 



Boltzmann equation are necessary to curtail the singu- 
larities at e = Ao + nui (n = 0,1,2...). However, so 
long as these singularities are localized on energy inter- 
vals that are much smaller than ui, the approximate Sn^ 
calculated from Eq. © will be suitable for the calcula- 
tion of both the enhanced order parameter via Eqs. [I] 
and [5] as well as the value of \ m the Ginzburg-Landau 
equation (Eq. |4} [8]. As expected from the analogy to 
Eliashberg's work, our singularities have energy widths 

of about Aj^^hp-jj. Hence, we may consider the in- 
equality — Qy*Q^!y./3 *C 1 as a further requirement for 
the validity of our linearization of the Boltzmann equa- 
tion for the calculation of \ m the superfluid case. We 
may also note that due to the fact that hui < 2Aq, no 
new quasiparticles are excited from the lower branch by 
pair breaking. Hence, the quasiparticle number is con- 
served in this first order approximation (X)k^ n k = 0). 
The quasiparticles are simply redistributed from the gap 
edge to higher energies. Substituting Sn^ into Eq. @, 
we find that at T ~ 0.13 Tp we calculate an increase in 
A by a factor of 1/10. So long as T < T c (0) , the relative 
enhancement increases with temperature because there 
are more particles to redistribute and the pulse does not 
have enough energy to break Cooper pairs. Unlike in 
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FIG. 3: The order parameter A as function of the temperature 
which solves the nonequilibrium gap equation for parameters f2(, 9 = 
110f2f , |q| = 0.1 kp, and hui = O.lbep. The black dashed line is 
the equilibrium dependence while the red dashed line gives the 



nonequilibrium transition temperature Tbcs > ^BCS' 
constrained Aq > hui 1 2 to avoid pair- breaking. 
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the normal fluid where Ao = at equilibrium, the en- 
hancement \ depends on the initial value of the gap. In 
Fig. |31 we plot the temperature dependence of the en- 
hanced nonequilibrium gap with a slightly stronger pulse 
given by f^g = 1100^, Tiuj — 0.15e^ , and |q| = 0.1 kp. 
Note that after the pulse is applied at equilibrium, the 
temperature of the bath may be increased above 
while maintaining a nonzero gap. This temperature de- 
pendence is exactly what would be expected from the 
analogous plot in Ref. Q. The new BCS transition tem- 
perature T c is given by the maximum of the nonequilib- 
rium plot where the inequality Tiui < 2Ao is saturated. 
If the temperature is further increased, A discontinu- 
ously vanishes again. For these parameters we calcu- 
late a small 3% increase of the transition temperature 
as expected from our requirement that nonlinear effects 
(even gap enhancing effects) be ignored in the Boltzmann 
equation. The temperatures we have considered are suf- 
ficiently close to the equilibrium transition temperature 

such that the approximation Tbcs ~ ^bcs(1 + x) is j us ~ 
tified. For small Ao a crude, order of magnitude estimate 

of x may be given by x ~ 2 ^""g r ° n ( A ) n ( A + hui). 



B. Normal at Equilibrium 

For contrast, let us now assume that we start in the 
normal phase at equilibrium. The initial distribution 
function is simply the Fermi-Dirac function for the nonin- 
teracting quasiparticles of Fermi Liquid theory. We have 
Ao = 0, so Eliashberg's requirement that tiui < 2Ao can- 
not be satisfied. However, we can guarantee that parti- 
cles are not excited from below the Fermi level by choos- 
ing ui and q such that the constraint 4e q £j? < (hui — e q ) 2 
is enforced. Because of this requirement, momentum and 
energy conservation cannot be simultaneously achieved 
for particles with energies less than ep. As such, only 
particles outside the Fermi sphere can undergo transi- 
tions. The lower equilibrium occupation number and 
higher density of states at high energies ensures that 
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FIG. 4: The T = 0.18 Tf normal phase occupation as a function 
of energy for Q, bg = 110 Qf, |q| = 0.1 fop, and Tiw = 0.21£f- 
The enhanced nonequilibrium transition temperature is Tbcs = 
1.3 rj^lg The dashed line is the occupation at equilibrium. 
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FIG. 5: The enhanced gap A as function of temperatures above 
^BCS ^ or f erm i° ns m the normal phase at equilibrium and param- 
eters f2;, 9 = 110 f!f, |q| = O.lkp, and hui = 0.21 sf 



quasiparticles just outside the Fermi sphere are excited 
to higher energies. 

Although in our system we have particle conservation 
just as in the superfluid case, we comment here that if 
these excited particles are allowed to leave the trap, then 
we have effectively cooled the fermions by sharpening the 
Fermi step, and the boson bath is unnecessary. Substi- 
tuting Sriu into Eq. we see that the depression of the 
population at the Fermi level shown in Fig. @] allows for a 

nonzero A above T B °^ S . As the temperature is increased, 
the nonequilibrium gap enhancement is overpowered by 
thermal smearing of the distribution function. There are 
more quasiparticles at energies near ef where they most 
strongly hinder superfluidity. This contrasts with the 
superfluid situation wherein the enhancement increases 



with temperature so long as T < T B °^ S . This effect can 
be seen in Fig.[5j where we find an enhancement of Tbcs 
by about 30% for parameters Slbg = HO Of, |q| = 0.1 hp, 
and huj = 0.2l£i?. This increase is much more dras- 
tic than the enhancement in the superfluid phase be- 
cause the requirements that Sn^ <C and huj <C £p 
are much less stringent than the superfluid requirements 
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— Q F ( 6 jv F )i/3 ^ 1 an( l ^ W — A- Thus, we may use a 
stronger pulse while still linearizing the Boltzmann equa- 
tion legitimately. As such, fermionic superfluidity is ex- 



pected to appear at temperatures as high as T : 
IV. SUMMARY 
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To conclude, we have shown that perturbing a system 
of trapped fermions creating a stationary quasiparticle 
distribution can be an effective way to stimulate fermion 
pairing and superfluidity. To demonstrate this, we cal- 
culate enhancements of the BCS order parameter and 
the transition temperature for a mixture of 87 Rb and 6 Li 
that is pushed out of equilibrium by a Bragg pulse. The 
mechanism by which fermions within the Fermi sphere 
are not excited differs depending on initial conditions. 
If the gas is a superfluid at equilibrium, these excita- 
tions are precluded by keeping Tito < 2Aq. In the normal 
phase, the parameters of the pulse can be chosen such 
that fermions below a certain energy cannot simultane- 
ously satisfy momentum and energy conservation. Thus, 
they are not excited. In both cases, the enhancements 
that we calculate are small, but this is a consequence of 
our perturbative treatment rather than a physical con- 
straint. This is evidenced by the strong effects reported 
from experiments on superconductors fl2l . liH ] , which are 
based on the same underlying mechanism. Finally, we 
suggest that by enhancing or creating a discontinuity in 
the quasistationary strongly non-equilibrium distribution 
of fermions (not necessarily at the Fermi momentum) via 
the technique proposed in this paper, one may achieve ef- 
fective BCS pairing at nominally very high temperatures 
of the bath. 
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